Joule Heating and Current-Induced Instabilities in Magnetic Nanocontacts 
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We consider the electrical current through a magnetic point contact in the limit of a strong inelastic 
scattering of electrons. In this limit local Joule heating of the contact region plays a decisive role in 
determining the transport properties of the point contact. We show that if an applied constant bias 
voltage exceeds a critical value, the stationary state of the system is unstable, and that periodic, 
non-harmonic oscillations in time of both the electrical current through the contact and the local 
temperature in the contact region develop spontaneously. Our estimations show that the necessary 
experimental conditions for observing such oscillations with characteristic frequencies in the range 
10* 10^ Hz can easily be met. We also show a possibility to manipulate upon the magnetization 
direction of a magnetic grain coupled through a point contact to a bulk ferromagnetic by exciting 
the above-mentioned thermal-electric oscillations. 
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INTRODUCTION: JOULE HEATING AND 
CURRENT-INDUCED INSTABILITIES IN 
MAGNETIC NANOCONTACTS. 



Electric transport in mesoscopic magnetic structures 
has become a hot topic of modern solid state physics re- 
search. One reason for this is the prospect for exploiting 
not only the charge of the electron but also its spin degree 
of freedom for electronics ( "spintronics" ) . 

An interplay between macroscopic magnetic degrees 
of freedom and the electron spin has already been dis- 
covered in the giant magnetoresistance (GMR) effect^'^. 
This is a spin-dependent tunnelling phenomenon that 
quickly found important practical applications in data 
storage devices. As one shrinks the dimensions of a mag- 
netic material toward the nanometer scale a new, to some 
extent inverse phenomenon occurs since a high density 
electric current is able to affect the magnetic order'^'"'. A 
current-induced precession and switching of magnetiza- 
tion has indeed been observed in a number of experiments 
on magnetic layers and point contacts^iSiL2i2ii£. 

Recentljiii a correlation was observed between the 
magnitude of a spin-torque effect and electronic scat- 
tering in a point contact between a normal and a fer- 
romagnetic metal. Using the method of point contact 
spectroscopy it was proven that elastic backscattering 
of electrons plays a crucial role in the transfer of mag- 
netization between differently magnetized regions of the 
ferromagnetic. This observation also raises the question 
of the role of inelastic scattering of electrons, including 
scattering by magnons, in magneto-transport since elas- 
tic scattering also shortens the inelastic diffusion length 
making the transport of heat away from the point contact 
less effective. 

An advantage of using electrical point contacts for 
switching the direction of magnetization is the possibility 
to locally achieve extremely high current concentrations 
of up to 9 10 A/cm^. Such high current densities are 



just what is needed for inducing magnetization dynam- 
ics. However, there is a restriction due to Joule heating 
that limits the possibility to increase the current density 
much further. This is because thermal heating eventu- 
ally destroys the contact. Nevertheless, for not too high 
temperatures, a special reversible point-contact trans- 
port regime arises in which electric and thermal transport 
processes leads to a spatial distribution of electric field 
and temperature that allows the temperature in the point 
contact to be controlled electricalljiiSiiiii. Heating may, 
however, significantly affect the resistivity of the material 
creating a nonlinear point-contact transport regime. As 
we will show below such a nonlinear response to an ap- 
plied DC voltage can lead to an electrical instability and 
to time dependent (oscillatory) transport phenomena. 

It is well known that strong Joule heating of a metal 
can result in an N- or S'-shaped current-voltage (J — V) 
characteristics containing segments with negative differ- 
ential conductance (see, e.g., the review Il5|) . In particu- 
lar. Joule heating of an anti-ferromagnetic break-junction 
was recently shown to result in an TV-shaped I — V 
curvei^. This phenomenon takes place if the balance be- 
tween the power released in the heated area and absorbed 
by the surrounding medium can be satisfied at three dif- 
ferent temperatures (two of them stable and the third 
unstable). As has been shown theoretically and experi- 
mentally, such a bi-stability in bulk metals at low temper- 
atures results in spontaneously formed thermo-electric 
domainsiiiSiiSiSSiSi, nonlinear periodic oscillations^^ and 
filaments^^ of high temperature and current in the sam- 
ple. This effect can be of prime importance for magnetic 
point contacts and particularly for those made of ferro- 
magnetic manganites. This is because the temperature 
dependence of their conductance is strong and diverse 
up to the Curie temperature and above, which may re- 
sult in N-, S'-shaped I — V characteristics as well as their 
combinations. 

In the present paper we investigate an electrical in- 
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stability caused by Joule heating of a magnetic point- 
contact junction. We show that non-harmonic periodic 
oscillations in time of the total current and the voltage 
drop across the Joule-heated part of the micro-contact 
(as well as the local temperature at the contact) appear 
spontaneously if the voltage bias of the entire sample is 
kept constant. This phenomenon takes place in a wide 
interval of sample temperatures that includes the Curie 
temperature of the ferromagnetic conductor. We also 
show that in the latter case the magnetization direction 
of a magnetic grain (which is coupled to a bulk ferromag- 
netic through a point-contact) periodically switches from 
the direction of the bulk ferromagnetic to the direction 
of the external magnetic field following periodic heating 
and cooling of the point-contact under the regime of the 
thermal-electric self-exciting oscillations. 



II. FORMULATION OF THE PROBLEM 

Below we consider the case when the electron relax- 
ation lengths for both momentum and energy are shorter 
than the size of the micro-contact (shown in Fig. 1). Un- 
der this condition the thermal and electrical characteris- 
tics obey the continuity equations for the energy 

rIT 

Cv^+d\v{-K{T)WT) = -a{T){W^)^ (1) 
and the local charge neutrality condition 



divj = 0, 



(2) 



while the total current / flowing through the system sat- 
isfies the equation 



Cdl/dt + RI = V 



(3) 



(see, e.g., Ref. |2^. Here Cy is the heat capacity of the 
metal per unit volume, T is the temperature, k is the 
thermal conductivity, a is the electrical conductivity, ip 
is the electric potential, j = —a{T)VLp is the electrical 
current density, C is the total inductance of the circuit, 
R is the total resistance (we assume below that the main 
contribution to R is from the point-contact) and V is the 
applied bias voltage. In writing Eq. we neglected the 
thermopower, which is small by a factor k^T/tp where 
fee is the Boltzmann constant and ep is the Fermi energy. 

The boundary conditions for the set of equations ^ 
and |2Jl are determined by the absence of heat and charge 
transport through the boundaries of the contact. It fol- 
lows that 

j„(r e E) = n-j = 0, jq,„(rel]) = n-jg = 0, (4) 

where jq = — KVT-|-i^j and n is the normal to the bound- 
ary E of the micro-contact. For a symmetric contact one 
has 




FIG. 1: Model of a micro-contact whose transport properties 
are discussed in the text. 



Here Tq is the temperature of the peripheral regions of 
the contact, the 2;-axis is directed along the axis of the 
contact. 

As was shown in Ref. it is mathematically conve- 
nient to use the coordinates u,v,(j) of an oblate ellip- 
soid for studying the kinetic properties of micro-contacts. 
These coordinates are related to the Cartesian coordi- 
nates X, y, z as 

X = do sin u cosh v cos , y ~ d^ sin u cosh v sin 
z — do cos u sinh v (6) 

where < u < 0, —oo < v < oo, < ((> < 2tt, u — 
6 = const at the boundary S of the micro-contact, c?o = 
d/(2sin0) is the effective length of the contact and d is 
the smallest diameter of the contact as indicated in Fig.^ 

We assume the contact to be symmetric and hence the 
temperature and the electric potential depend only on 
the coordinate v and time t. Therefore, Eqs. 1^ and |(5J, 
expressed in the new coordinates, reduce to 
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^^dT 

dt dl sinh^ v cosh v dv 



^0 

a{T) 



dT 

k(T) coshw— — 
ov 



d_ 

dv 



d^ sinh^ 



cr(T) cosh 



dip 
V \ dv 
d(p 
dv 



= 



(7) 
(8) 



While deriving Eq. © we neglected the term cos^ u in 
the sum sinh^ v + cos^ u (which arises in the Lame coeffi- 
cients). This is a reasonable approximation if the charac- 
teristic size of the region of Joule heating is large enough. 

Integrating Eqs. lO, lO and the relation j = —a{T)Vip 
over the cross-section area of the sample one gets a set 
of equations that describe the space- and time-evolution 
of the temperature T{ri, t) and the total current I{t) of a 
Joule- heated microcontact, 



p{T)I^/ll = 0; (9) 



y 

Lp{z ±oo) = ± — ; 



T{z ±oo) 



(5) 



i?=(p(T)) 



(10) 
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In Eq. where f{i]) = dp sinh^ 2u(ry)/4 and = 
27r(iQ(l — COS0), we introduced a new variable 77 that is 
related to u through the relation 



d , d 

— = coshu— , 
or] ov 



(11) 



while in Eq. (|10l) and below, angular brackets imply an 
integration over 77, 
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-7r/2 



^0 J-TT 



t/2 



.drj 



III. NEGATIVE DIFFERENTIAL RESISTANCE 
AND STABILITY OF THE STATIONARY 
TEMPERATURE AND CURRENT IN A 
JOULE-HEATED MICROCONTACT. 

The set of equations JHl and 1)10(1 subject to the bound- 
ary conditions Q always have steady-state solutions 
/ = i{V) and T = T{V, 77) that satisfy the stationary 
equations 



d_ 

d-q 



K{f) 



df 
drj 



p{f)P/ll = 



V = RI 



(12) 



(13) 



According to Refs. |2a and |23 the stationary current- 
voltage characteristics of the microcontact is 



i{v) 



V 



R{vy 



R- 



'(^) = 7 / ' dx/p (t,^{V)VT^) 



(14) 

where Tm is the temperature in the center of the mi- 
crocontact (77 — 0) and depends on the applied voltage 
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rr\2, /T~i2 

-'m — -'bulk 



4:L 



(15) 



where L = 7r^fc^/(3e^) is the Lorentz number. The de- 
pendence of the stationary temperature T{V, 77) on the 
bias voltage V and the coordinate rj is determined by the 
relation 

RiV) r^^^ dx/p {tUvWi~x^) - ^ (16) 



where A(T) = ^ {AL/V^){T^ - T^^,^) 

Differentiating Eq. ((14(1 with respect to V one sees 
that the differential conductance dl/dV is negative if 
-R{V){VdR~^/dV) > 1. Assuming the resistivity to 
be p{T) = po -t-pi(r), where po is the residual resistivity 
and the temperature-dependent part pi oc T", one finds 
the differential conductance to be negative if 



1 - R 



Po 



dx 



p{T,^VT^)J p{T„,VT^) 



> 1 



(17) 



If the microcontact is heated to temperatures for which 
Pi{Tm) ^ Po the inequality ifTTI) reduces to 



a > 1 



Po 



(18) 



where fc is a constant of order one (k ^ 1). For many fer- 
romagnetic metals a ~ 2 and increases as the tempera- 
ture approaches the Curie temperature (see, e.g., Ref. l29j) 
and hence the I — V characteristics of a Joule-heated 
ferromagnetic microcontact has sections with a negative 
differential resistivity in a wide range of temperatures 
including the Curie temperatures. 

As soon as the micro-contact is heated to temperatures 
where the differential resistance is negative, the station- 
ary distribution of temperature and current in the con- 
tact may be unstable. 

In order to investigate the stability of the stationary 
temperature T and current / given by Eqs. (EH and lO 
we write the temperature T(i, 77) and the current I{t) as 
sums of two terms, 

r(i,7;) = r(77) + ri(t,7;), /(t) = / + /i(t) , (19) 

where Ti and Ii are small corrections to the stationary 
values. After substituting Eq. JTHl) into Eqs. © and jllHl 
we obtain a linearized set of equations for Oi = k{T)Ti 
and Ii of the form 



dt 



lo 



(20) 



/:§ + i?/, + .i?/ = o; ^i?-(^e,^, (21) 

where /? = l/(/(?7)Cy(T')) and the Hcrmitian operator is 



H 



_ I^P{T) 
drj^ ~ lln{T) 



(22) 



(here and subsequently, a prime denotes differentiation 
with respect to T). 

As shown in Appendix fXl the Laplace transformed cur- 
rent 

/>oo 

ii{p) = i I{t) exp{—pt}dt 
Jo 



can be expressed as 



7i(p) = {/i(0) - i/^ A^}/Dip) , (23) 



where the denominator is 




C ll ^ K{v + K) ] Cdl 
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Here are coefficients of the series expansion of 81 
in eigenfunctions V'l/ of the Hermitian operator Hp — 
(3{v)H, and \,j are the eigenvalues of this operator (see 
Eq. E^ ). 

It follows from Eqs. H23|l and H24|) that for the station- 
ary current / and temperature T (see Eqs. (|14ll and (|16() ) 
to be stable it is necessary and sufficient that the func- 
tion D{p) has no zeros in the half-plane Re p > 0. It is 
easy to see that the stationary current and temperature 
are stable if the inductance C is small, that is tj <C tt 
where 



TT 



P p{T) 



(25) 



are the characteristic relaxation times for current and 
temperature, respectively. In this case D{p) = at p cx 
1/t/ and hence one can expand Dij)) as a power series in 
and find Ii{t) cx exp {pi^2t) where 



Pi, 2 



(26) 



that is the stationary distribution is stable at low induc- 
tances. 

In the opposite case of large C (r/ ^ r^) one can 
neglect the second term in Eq. 124|l and find 



A. Adiabatic evolution equations 

In order to investigate the evolution of the instabil- 
ity (which is governed by Eqs. and we assume 
that the inductance C is large enough to have tj ^ tt- 
In this case the current varies slowly in time while the 
temperature rapidly follows the current varations. This 
assumption allows us to develop an adiabatic perturba- 
tion theory in which the temperature and the current can 
be represented as 



m = i{t) + h{t). 



(28) 



Here T{U{t),j]) is the steady-state solution of 
Eqs. 112(1 and (|13() in which the constant voltage V has 
been replaced by a time-dependent voltage U{t) (to be 
determined). Hence T{U{t),ri) satisfies the equation 



d_ 

drj 



~ dT 



^p{f)P{U)/ll^O 

I{U) = U/R{U) , (29) 



where 



R{U) = R{f} = (pit) 



(30) 



, IdV 

Ii{t) (x exp \ 



(27) 



From Eq. (|27|l it follows that the stationary current and 
the temperature are always stable in the branches of the 
I ~ V curve that have a positive differential resistance, 
and they are unstable in the branches with a negative 
differential resistance. In the next section we investigate 
the adiabatic evolution of this instability in the case of a 
large inductance C. 



while the current I{t) satisfies the equation 



dl ~ ~ 

C— + R{T}I = V- 



(31) 



An equation for U (t) can be found from the condition 
that Ti and Ii arc small corrections to T and /, respec- 
tively (that is iTil < f and < /). As shown in 
Appendix IbI this condition is satisfied if the current I{t) 
and the voltage U{t) satisfy the following set of ordinary 
differential equations: 



IV. THERMO-ELECTRIC SELF-EXCITED 
PERIODIC OSCILLATIONS OF THE TOTAL 

CURRENT AND OF THE VOLTAGE DROP 
OVER THE JOULE-HEATED PART OF THE 
MICROCONTACT 



As was shown in Section iHIl the stationary current / 
and temperature T{ri) in a microcontact (Eq. H14|) and 
Eq. become unstable if r/ > tt- Below we show 

that this instability results in the spontaneous appear- 
ance of non-harmonic periodic oscillations in time of the 
total current / flowing through the microcontact and of 
the voltage drop U over the Joule-heated part of the con- 
tact provided the applied bias voltage V is kept fixed. 
In Subsection IIV Al we derive a set of reduced adiabatic 
equations for / and U; in Subsection IIV Bl we study so- 
lutions of these equations both analytically and numeri- 
cally. 



C^+R{U)i^V; 



(32) 



Here 



,=4i^?)-'r^^) ™ 



dU J \ dU k 



R^ 



and V is the applied bias voltage; the I — U curve 
I{U) = U/R{U) and the resistance R{U) are determined 
explicitly by Eq. H14|l and Eq. (|15l) in which V has to be 
changed to U{t). 

The set of non-linear ordinary differential equations 
Ea. ((32|l allows the system under consideration to be 
mapped onto the effective circuit shown in Fig.(|2Jl. The 
circuit includes a Joule-heated conductor in series with 
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FIG. 2: An effective circuit under a fixed voltage V containing 
a Joule-heated conductor in series with an inductance C; the 
resistance R{U) depends on the voltage drop U on the con- 
ductor in such a way that the conductor CVC I{U) = U / R{U) 
has a section with a negative differential resistance; I is the 
total current flowing in the circuit. 



an inductance £, the total voltage drop in the circuit 
V being kept constant in time. The Joule-heating of 
the conductor is assumed to be inhomogeneous, so that 
the relationship between its temperature and voltage U 
is described by the thermal continuity equation Ea. H28|) 
and hence the rate of heat removal from the conductor is 
Q = RP{U). 

In the next Subsection we investigate the conditions 
under which the steady-state solution of Eq. (|32() is un- 
stable and how it evolves into a limit cycle in the I — U 
plane, a limit cycle that corresponds to periodic non- 
harmonic oscillations of the current and the voltage drop 
over the Joule- heated part of the micro-contact. 



B. Spontaneous development of periodic 
oscillations of the current and of the voltage drop 
over the point contact 

Let the steady-state solution of the set of equations 
(|22J) be 



Jo = /(f/o); Uo = V 



(34) 



A study of the stability of the steady-state solution 
with respect to small perturbations, carried out on the 
basis of the linearized set of equations (|32|l shows that 
small variations Ji and Ui from the steady-state develop 
as exp (71. 2^) with time, where the rate factor 



71,2 



1 dl I 
C^^dVji 



± 



'1 dl I 

^c^^dvjl 



21 

(35) 

From Eq. I|35|l one sees that in the case of a negative 
differential conductance the steady-state solution looses 
its stability with an increase of the inductance as soon as 
the term in the round brackets changes sign, i.e. as soon 
as £ > Crr, where the critical value of the inductance is 




tU/Uc 



FIG. 3: Dependence of the total current / through the micro- 
contact shown in Fig. 1 on the voltage drop U over the con- 
tact for the case that the micro-contact resistance depends on 
voltage as R{U) = {Ro + PU^). The current is normalized to 
-^0 = Uq/Ro, where Uo = \/Ro/P 



Comparing Eq. H36|l with Eq. H25|) one sees that sta- 
bility is lost as soon as r/ > tt in accordance with the 
above analysis. Therefore, if tj > tt any initial state 
close to the steady state Uq, Iq is repelled from this point 
in the voltage-current plane. On the other hand, as one 
easily sees from Eq. 1321 any initial state which is very 
far from the stable point (that is |/(< ~ 0)| ^ |/o| or 
\U{t = 0)1 3> \Uo\) decreases in time and is attracted to 
the stable point. It means that there is a stable limit 
circle in the I —U plane (see, e.g., Ref. IIO)- This means 
that non-linear periodic oscillations of the total current 
and the voltage drop over the micro-contact appear spon- 
taneously if the bias voltage V is kept fixed. In order to 
show characteristic features of the spontaneous develop- 
ment of electric oscillations, we present a typical current- 
voltage characteristics for a micro-contact in Fig. 13 and 
the corresponding spontaneous oscillation cycle of the 
current through the system and the voltage drop over 
the Joule-heated micro-contact is shown in Fig0| 

Simple estimations of the ratio tj/tt (see Ea. (|25f) ') 
show that the critical value of the inductance Ccr (at 
which ti/tt ~ 1 ) is 



TCy 

Jo do 



(37) 



where jo is the characteristic value of the current den- 
sity in the micro-contact and do is the size of the micro- 
contact. 



C. Magnetization switching under the 
thermal-electric self-exciting oscillations. 



R 



dV_ 
dl 



(36) 



The system under consideration is presented in Fig|Sl 
A conducting magnetic grain is coupled through a point- 
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FIG. 4: Spontaneous oscillations of the current I{t) trough 
the micro-contact sketched in Fig. 1 and the voltage drop 
U{t) over the contact. Here Jo ~ Uo/Ro and Uo = \/Ro/f3. 
Initial values of /(O) and ?7(0) are shown by filled squares. 
As time goes on the current and the voltage drop move along 
the dashed or thin solid line towards the limit cycle (thick 
solid line) depending on whether the initial state is outside or 
inside the limit cycle. After a time t 3> t/ the set {I{t), U{t)) 
moves along the limit cycle executing a periodic motion, i.e. 
the current I{t) and the voltage drop U{t) execute periodic 
non-linear oscillations. 



contact (PC) to a bulk magnetic conductor and to a non- 
magnetic conductor (the right-hand side of the figure). 
There is a magnetic field H directed opposite to the mag- 
netization of the bulk magnetic conductor. The magnetic 
field is weak enough so that at low temperatures the grain 
magnetization is parallel to the magnetization of the bulk 
ferromagnetic due to the exchange interaction. 

Let the inductance of the circuit to be large enough 
in order that the above-mentioned thermal-electric self- 
oscillations arise in the system. The limiting cycle along 
which the current J and the voltage drop U periodically 
move in time are shown in Fig|Sl(for the sake of simplicity 
it is shown for the case of a extremely large inductance 
T£ > tt). 

When J and U move along the left arising branch of 
the CVC the PC is cold (its temperature is nearly the 
same as those of the cooling media). When the system 
jumps to the right arising branch of the CVC (along the 
upper straight line in FigEJ, the PC is Joule heated above 
the Curie temperature while the bulk ferromagnetic and 
the grain remain cold. It means that the magnetization 
in the vicinity of the PC disappears and hence the ex- 
change coupling of the grain to the bulk ferromagnetic 
is interrupted, and the magnetization of the grain fiips 
to the direction of the magnetic field H. Moving further 
along the limiting cycle (along the right branch of CVC), 
the system jumps to the left arising branch of the CVC 
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FIG. 5: A. The nanocontact is cold. The exchange coupling 
through the contact keeps the magnetization of the grain par- 
allel to the magnetization of the left bulk ferromagnetic con- 
ductor. B The nanocontact is hot. The exchange coupling 
through the nanocontact is killed and the magnetic field H 
fiips the magnetization of the grain. 



Uo 



U 



FIG. 6: The limiting cycle for an extremely large inductance 
in the circuit {tc S> tt); when the total current I{t) and the 
voltage drop U{t) on the point-contact move along the left 
branch of CVC the point contact is cooled to temperatures 
close to the thermostat temperature; while they move along 
the right branch of CVC the point contact is heated above 
the Curie temperature. 



(along the lower straight line in Fig|SJ), the PC is cooled 
to temperatures close to the thermostat temperature, the 
magnetization in the vicinity of the PC is restored, and 
the grain magnetization fiips to the direction of the bulk 
magnetization again. After that these fiips repeat period- 
ically in time with the periodicity of the thermal-electric 
self-exciting oscillations. 

The period of the thermal-electric self-oscillations and 
hence the period of the switching of the grain magneti- 
zation is controlled by the inductance in the circuit (as 
well as by the applied voltage). 
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The above-mentioned mechanism of the magnetization 
switching works if, at low temperatures, the exchange 
coupling of the grain to the bulk ferromagnetic keeps the 
grain magnetization parallel to the bulk magnetization 
despite the opposite direction of the applied magnetic 
field. By this is meant that if at low temperatures, the 
grain magnetization flips to the direction of the magnetic 
field the total energy of the system increases, that is 

Wtotai = I dV\a (^^^ - Mi/j > (38) 

Here M is the magnetic moment, a ~ I/aM^ where 
/ ~ kTc and a are the exchange energy and the lattice 
spacing, respectively, k is the Boltzmann constant and 
Tc is the Curie temperature. If the characteristic size of 
the PC d is less than the domain wall length Idw, the 
length of the effective domain wall is of the order of d 
(see^i), and from Eq.(|2HIl it follows that Wtotai > if 
a{M^/d'^)d^ - {hiOH/a^Wo > that is 



Vo<^ (a'd) 



(39) 



where Vq is the volume of the grain, Th — hLUn/k and 
luh = eH/{mc) (e is the electron charge and m is the 
electron mass). If one takes, e.g., d ~ 10a, Tc ^ lO'^K, 
H - 10^ lO^gauss (that is IQ-^ IQ-^ T)and Vb 
10a X one finds the characteristic width of the grain 
L ~ (10^ ~ 10'^)a where a is the atomic spacing in the 
crystal. 



V. CONCLUSION 

In conclusion we have theoretically studied the tem- 
perature dependent electrical current through a voltage- 
biased point contact. We predict that spontaneous tem- 
poral oscillations of the temperature and current through 
the contact may appear due to Joule heating and a strong 
temperature dependence of the point contact resistance. 

We have also showed that the self-exciting oscillations 
can control the magnetization direction of a magnetic 
grain coupled through a point-contact to a bulk mag- 
netic conductor: the grain magnetization direction pe- 
riodically switches from the direction of the bulk mag- 
netization to the direction of the external magnetic field 
following periodic cooling and heating of the point con- 
tact under the regime of the self-exciting thermal-electric 
oscillations. In order to estimate the characteristic fre- 
quency of these self-oscillations we take the minimum 
width d of the micro-contact shown in Fig. 1 to be 200 nm 
and let the current / through the contact be^^ 1 mA, its 
temperature T '--^ 30 K and the specific heat Cy at this 
temperaturaSSi ~ 1 Jcm~'^K~^. Under these conditions 
we can use Eq.l^Hlto evaluate the critical inductance £cr 
required for self-oscillations to develop as well as the os- 
cillation frequency lu when C > Ccr- The result is 



and we note that the frequency of the self-oscillations 
decreases linearly with an increase of the inductance of 
the circuit. 

Finally we observe that due to Joule heating of the 
contact area, where the current density is maximal, a 
wide range of temperatures can be reached in this region 
of the sample (and controlled by the bias voltage). This 
allows for a thermal scanning that avoids difficult calori- 
metric measurements in many cases and since the Joule 
heating easily can produce temperatures of the order of 
the critical temperature for magnetic phase transitions, 
it offers an exciting possibility for electrically controlled 
magnetic switching on the nanometer scale. 
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APPENDIX A 

Let us write 61(77, t) in the form of an expansion as 
ei(7/,i)=5]A.(t)V.(ry), (Al) 

V 

where tl}^ are the eigenfunctions of the operator Hp = 
(3{ri)H, which satisfies the Sturm-Liouville equation 



lln{f) 



(A2) 



Substituting Eq. (|Xl1) into Eqs. lP7|) - (|B3)l and carry- 
ing out a Laplace transformation with respect to t we 
find the following set of algebraic equations for a(p) = 
A{t) exp {—pt)dt and ii{p) = I[t) exp {—pt)dt: 



{p + K)a,-{2I/ll)p,ii = A,(0) 
{Cp + R)ii+I^R^,a^ = Ch{Q). 



(A3) 



Here Ai, (0) and Ii (0) are the initial values of A^, and Ii 
at t = 0, respectively, while 



R,, — 



p{T) 



^•^iv) ) ; Pi^ = {p{t)iIjAv)) 



While obtaining Eq. (|A3(I we took into account the fact 
that the Hermitian conjugate operator i?^ = HP{ri) has 
eigenfunctions ipi^ — f3~^ip,y and hence {ip^(l(3)ip^>) — 

On the other hand, differentiating the current-voltage 
characteristics V = (^p{T{I ^rf)^ I and Eq. (|12|l with re- 
spect to the current / one obtains 



and 



10-^ 



(40) 



dV „ t/ 'dT 



(A4) 



(A5) 
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Here H is the operator defined by Eq. H22() . It follows 
from Eq. (|A5|I and Eq. (|A2() that ndf/dl can be written 
in the form 



df 2/ ^ p. 



Substituting Eq. ijMjl into Eq. (|X4|) we find 



dl 



if) 



(A6) 



(A7) 



Using Eq. IjASp and the equality Eq. (|A7p one gets 
Eq.(|^ 



APPENDIX B 



Inserting first Eq. into Eqs. © and l(Tn|) we obtain 
linear differential equations for Ti and /i of the form 

Hf,(i^fi)-2/3pif)^h = PF (Bl) 



where 



C^+R{f}h+SRi = 0. 
dt 



6R = (p'{f)f,)- P^p{f). 



(B2) 



(B3) 



When deriving the above equations we neglected a term 
oc dTi/dt in Eq. ljBl|l since it is of second order in the 
parameter 7 = tt/tj <C 1 (in contrast to Cdli/dt, which 
is of the first order in 7). 

Equation (|B1|I can be easily solved if one writes 9i = 
k(T)Ti in the form of a series expansion as 



(B4) 



Using now Eq. (|A7|) (in which all quantities marked with 
the bar sign should be changed to those marked with the 
tilde sign) one finds the equation for h to be 



dh , dU~ ^R,F, 



A. 



(B7) 



Solving the linear differential equation ljB7|) one finds the 
solution of the linear differential equations ljBl|l and ljB2p 
for the current Ii{t) to be 



h{t) ^ exp{r{t)} 



cxp{-r(T)}+/i(0)} (B8) 



where 



(B9) 



and /i(0) is the initial value of /i(i). From Eqs. ljB8|) and 
(|Ij9|) it follows that the current Ii {t) grows exponentially 
with time if the differential resistance is negative (i.e. if 
dU /dl < 0). In order to prevent this increase we impose 
a condition on the parameters of the zero approximation 
incorporated in F^. The condition is expressed by the 
equation 

-ZT'^^T.^ ^"^^'^ + (0) - (BIO) 
As a result one gets 

^ 1 em r dr V ^ e-^M (BU) 

^ Jt ,, \v 



It follows from Eq. (|B10|) and Eq. IjBlip that the inequal- 
ity <C |/| is satisfied at all times if 



Inserting the solution of Eq. (|Bip found in this way into 
Eq. HB2I) one may rewrite this equation to read 

Here '01/ and Ai, are eigenf unctions and eigenvalues of the 
operator (see Eq. HA2|l ) in which T and /(F) are 
changed to T and /(t/), while 



A. 



Therefore, the parameters of the zero approximation /(t) 
and Uit) should satisfy the equation 



(B12) 



F, = ( F{t, r;)0. ); Ru = / W = ( /^(^)V'- 



(B6) 



Using Eq. ||bT2)| and Eqs. JTHl, (jSzl, EU, and ^ we 

find the final set of equations Eq.(j221l for the current I{t) 
and the voltage U{t). 
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